THE ENERGY OF CRUMPLED SHEETS IN FOPPL-VON 
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Abstract. We study investigate a long, thin rectangular elastic membrane that is 
bent through an angle 2a, using the Foppl-von Karman ansatz in a geometrically 
linear setting. We study the associated variational problem, and show the existence 
of a minimizer for the elastic energy. We also prove rigorous upper and lower bounds 
for the minimum energy of this configuration in terms of the plate thickness a and 
the bending angle, and we also obtain results for the structure of the elastic ridge 
along it's length. 



1. Introduction 

Crumpled elastic sheets can be thought of as minimum energy configurations for 
the Foppl - von Karman (FvK) energy. Using this approach of elastic energy mini- 
mization, the crumpling response is understood as a result of the elastic energy of the 
sheet concentrating on a small subset of the entire sheet [21 123 1211 ■ The energy in 
a crumpled sheet is concentrated on a network of thin line-like creases (ridges) that 
meet in point-like vertices. Recent work has resulted in quantitative understanding 
of both the vertices □ |H1 El EH and the ridges [2HI 123 IHl US • Scaling laws gov- 
erning the behavior of crumpled sheets have been obtained j21 123 I2H| using scaling 
arguments. 

Minimum energy configurations for the FvK energy have also been studied in the 
context of the blistering problem, viz. the buckling of membranes as a result of 
isotropic compression along the boundary [221 HI] • 

There is a considerable body of mathematical work focused on the blistering prob- 
lem 1^ 123 13 1221 12ni El HI- Upper and lower bounds have been obtained for ap- 
proximations to the elastic energy |321 1^ E^i , for the FvK energy [221 and for full 
three dimensional nonlinear elasticity |3] . The FvK energy and full three dimensional 
nonlinear elasticity give the same scaling for the upper and the lower bounds. 

Our goal in this paper is to prove corresponding rigorous results for the energy in 
a minimal ridge - a single crease in a crumpled sheet. In addition to scaling results 
for the energy, we also investigate the structure of the ridge by obtaining pointwise 
bounds for its "width". Our results for the ridge show an interesting contrast with 
the corresponding results for the blistering problem |21 11] • In particular, the scaling 
of the energy with the thickness of the sheet has a different exponent. This implies 
that the boundary conditions play an important role in determining the F-limit of 
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the FvK energy in the hmit the thickness goes to zero. We discuss this issue further 
in Sec. 13 

This paper is organized as follows - In Sec. we describe the problem of interest, 
set up the relevant energy functional and determine the appropriate boundary con- 
ditions. We also rescale the various quantities to a form that is suitable for further 
analysis, and recast the problem in terms of the rescaled quantities. In Sec. HI we 
prove a lower bound for the elastic energy for our boundary conditions. In Sec. 
we prove the corresponding upper bounds by explicit construction of a test solution. 
In Sec. we investigate the structure of a single ridge, and we present a concluding 
discussion in Sec. [7| 

2. The variational problem 

We are interested in a minimal ridge, i.e., the single crease that is formed when 
a long rectangular elastic strip is bent through an angle by clamping the lateral 
boundaries to a bent frame. This situation is depicted in Figure [H 




Figure 1. A minimal ridge. The boundary conditions are given by 
a frame (the thick solid lines) bent through an angle. The sheet is 
essentially flat outside the region bounded by the two dashed curves, 
and the bulk of the energy is concentrated in this region. 

As we will see below, the idealized boundary conditions with a sharp corner are not 
appropriate, since they lead to an infinite energy for a sheet with a finite thickness. 
If we make the corner extremely sharp, all the energy (asymptotically) will be at the 
corners, and this obscures the interesting physics in the problem, namely the energy 
and the structure of the ridge. Thus we have to incorporate the smoothness of the 
corners into our boundary conditions. In general, curvatures on scale smaller than the 
thickness cannot occur for a real sheet, and our model energy is not appropriate for 
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this situation. For a crease in a real crumpled sheet, the corner is definitely smooth 
on the scale of the thickness of the sheet. 

We will consider the situation where the sheet is clamped to a frame, that is 
much like the idealized situation depicted in Fig. ^ The sheet is a rectangular strip 
\x\ < L,\y\ < L' . We generally consider the situation L < L'. The sheet is clamped 
to the frame at x = ±L and a.t y = ±L'. We will place the following requirements 
on the frame - 

(1) The boundary conditions at the frame are non-stretching, i.e, the strain 'jyy 
is identically zero on the boundaries x = ±L, and the strains j^x = Ixy = 
at 1/ = ±L'. 

(2) The bending at the boundaries is localized. If s is an arclength parameter 
along the boundary, we will require that both the boundaries are straight on 
a set of the form |s| > k, where k <^ L. 

(3) We will assume that the two "bent" boundaries on the frame {x = ±L) are 
given by planar curves, and that the angle between the straight sections on 
these boundaries ioi y > k and y < —k is 2a for both boundaries. 

(4) We will require that the straight segments for both the "bent" boundaries are 
parallel and a distance L apart. This implies that the two boundaries are not 
twisted with respect to one another. 

(5) We will also require that the two boundaries have the same asymptotic shift. 
This is necessary to make 'y^y = at y = ±L'. To define this precisely, we 
need to introduce appropriate coordinates, and we will do this below. 

2.1. Coordinates. We will use (the material) coordinates {x,y) on the reference 
half strip |x| < L, \y\ < L' . The planes containing the parallel straight portions of 
the two boundaries pick out preferred in-plane directions for y ^ k and y <^ —k. 
We introduce two sets of Cartesian coordinate systems in the ambient space. The 
coordinate direction u is perpendicular to the straight portion of the boundary at 
X = —L , and is directed toward the boundary at x = L. The coordinate directions 
f + and v~ are along the straight portions of the boundaries for s > k and s < —k 
respectively. The coordinate directions and w~ give the out of plane directions, 
and are chosen so that {u, f +, w^) and {u, v~ ,w~) are right handed orthogonal triads. 
Finally, the origins of both the coordinate systems coincide, and they are chosen 
such that the straight portions of the boundaries lie in the plane = 0, and the 
frame boundaries are in the planes u = ±L. Note that this completely specifies the 
definition of the coordinate systems, and in particular, we do not have a freedom to 
translate v^. 

The various coordinate are represented schematically in Figure El The grid in the 
figure is generated by the lines x = constant and y = constant that are straight in 
the reference (material) coordinates. We will use the coordinate system {u,v~^ ,w^) 
for the portion of the sheet with y >0 and {u, v~ ,w~) for y < 0. At |/ = 0, we have 
the matching conditions 

f w~ \ _ f cos 2a; —sin 2a \ / \ 
^ ^ \^ f~ y \^ sin2a cos 2a J \ ^'^ J 
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Also, in the straight portion of the boundaries |s| > k, since = 0, it follows 
that v'^{±L,y) — y is a. constant for sufficiently large (small) y. For the boundary at 
X — —L, we define the asymptotic shifts by 

Sf — v^{—L, y) — y, y > k (respectively y < —k). 

Similarly, the asymptotic shifts for the boundary at a; = L are constant if \y\ is 
sufficiently large. For the boundary at x — —L, we define the asymptotic shifts by 

6f — v'^{L, y) — y, y ^ k (respectively y < —k). 

We will say that the two boundaries are compatible, if 5f — 5j and 5f = 5^. This 
clearly a necessary condition for the existence of a configuration of the sheet that 
satisfies the boundary conditions, and is asymptotically strain free, i.e, the strain is 
identically zero for \y\ > / for a sufficiently large /. In particular, we can take / = k. 

Assuming that the two boundaries are compatible, we will set S'^ — — S2, 
S~ — — S2, and S — S'^ — S~. The quantities S'^ and S~ change under translations 
of the coordinate y, but S is an invariant under these translations, and is purely a 
geometrical property of the frame. Below, we will give an expression for S in terms 
of functions specifying the boundary conditions. W.L.O.G, we can, and henceforth 
will, translate the coordinate y such that 5+ = —6~ — ^S. 




Figure 2. A schematic representation of our coordinate system and 
the boundary conditions imposed on the sheet. The thick lines depict 
the "frame", and the "corners" are smooth on a scale a. and 
are the in-plane directions and are the out-of-plane directions. The 
grid is given by the lines x — constant and y — constant 



2.2. The Elastic energy. A mathematically justified way to obtain the elastic en- 
ergy of the deformed sheet is to treat the sheet as a three dimensional (albeit thin) 
object and use a full nonlinear three dimensional elastic energy functional for the en- 
ergy density. The sheet is now a three dimensional object S x [— |, |] with thickness 
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h. Let X, y denote the in plane coordinates as above, and z denote the coordi- 
nate in the thin direction. If the configuration of the sheet is given by a mapping 
(j) : S X [— |, |] — ^ the elastic energy is given by 

f f rh/2 

I3D= // dxdy / dzW3D{^<f)) 



JJs J 

This approach however does not take advantage of the "thinness" of the sheet. 
In particular, we would like to treat the thin sheet as a two dimensional object. 
The derivation of reduced dimensional descriptions of thin sheets has a long history. 
There is a classical theory for thin elastic sheets built on the work of Euler, Cauchy, 
Kirchoff, Foppl and Von Karman |3m |2Sl IHj • 

In the classical Foppl - von Karman ansatz, the behavior of the deformation is 
completely determined by the behavior of the center-plane z = 0. An asymptotic 
expansion with this ansatz P] yields an effective 2-D elastic energy 



X=h 



jj dxdyW2D{(t)x,4>y) + qh^ j j dxdy\Vv\ 



where g is a nondimensional 0(1) factor, v = ipx ^ 4'y/\\4>x x 4'y\\ is the normal 
to the center surface, and W2D is an effective two dimensional energy. This is the 
geometrically nonlinear Foppl- von Karman energy of the thin sheet. 

The functional W2D is zero if {(f)x,(j)y) G 0(2,3), the set of matrices that give 
isometric linear mappings — M'^. We will demand that 

W2D{<l)x,<Py) > cdist2((0,,,0^),O(2,3)). 

and that W2d{(Px, <Py) ^ ||0aH^ + ioi large 4>x, 4>y These conditions are identical 
to the conditions on the energy in Ref. 

A typical (or canonical) energy functional which satisfies these conditions, and has 
the natural invariances for the problem, viz., the action of 0(3) on {u, v, w) and 0(2) 
on (x, y) is 

iy2D = dist2((0,,0j,),O(2,3)) 

= {ul + vl + wl- 1)^ + 2{uxVx + UyVy + w^Wyf + {u^ + Vy + - 1)^ 

For the most part, we will work with the energy that we get by linearizing the 
above expression, only in the in plane deformation {u, v) about the reference state 
Ux = l^Vy = 1 and all the other derivatives are zero. This yields the linearized energy 

W2D = [2{ux -l)+wlY+ [2(^4/ - 1) + wij] V 2 [ux + Vy + WxWyf . 

Note that, by linearizing the energy functional, we have destroyed the natural in- 
variances of the energy, and have picked out preferred in-plane and out of plane 
directions in the ambient space We will also linearize the normal vector so that 
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Using this expression for the normal vector, gives the geometrically linear Foppl - 
von Karman energy 



s 



+2 [u^ + Vy + w^Wyf + qh"^ [wl^ + 2wly + wJJ ] 

We will henceforth normalize the energy by 4/i, and also define the scaled thickness 
a = \f\h. The linearized FvK energy for our problem is X = X"*" + X~, where 



j/=0 J x=—L 



1 1 

UX + - 1 j + - (f X + Uy + W^Wyf 



1 \ (J 



dxdy, 



where we have suppressed the subscript + on and w^. I~ is given by the corre- 
sponding expression for y < 0. 

As defined above x and y are reference coordinates on the sheet, u and v are in- 
plane coordinates, w is the out of plane displacement and a is the scaled thickness 
of the sheet. The integrand includes the squares of the linearized strains, 

IXX = - 1, 'Jxy = 7yX = ^{Vx + Uy + W^Wy), 'Jyy =Vy + " 1. 

2.3. Boundary and matching conditions. The blistering of thin films is also 
described by the elastic energy in (j21). A similar energy also describes multiple scale 
buckling in free elastic sheets {i.e. sheets that are not forced through the boundary 
conditions) that are not intrinsically flat I34j . 

The difference between the blistering problem and a minimal ridge is in the bound- 
ary conditions, which we describe below. If the bending half-angle a <^ 1, tana ^ 
sino; ~ a. In this case, the deformations and the linearized strains are small. 

If s is an arclength parameter along the boundary, with our choice of coordinate 
systems, the frame is given by functions {vf{s),wf{s)) for s > (resp. s < 0) at 
X = — L, and functions ivf{s),wf{s)) for s > (resp. s < 0) at a; = L. We define 
the curvature ni of the boundary at x = —L by 
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Ki = — = / {w'^ssYds + {w ssfds 

O'l Jo J~L' 



ai is a length scale associated with the boundary at a; = —L. We define ki by 
demanding that the boundary at a; = —L is straight for \s\ > ki. We will require 
that ki be on the same scale as ai, so that ki = Kai, with K staying order 1 as we 
change the parameters in the problem. We can similarly define 02 and K2 using the 
boundary at x = L. We will assume that the length scales Oi and 02 are comparable, 
i.e 02/01 ~ 0(1) and likewise for ki and /c2- In fact, we will typically assume that 
they are equal ai = 02 = a, and ki = k2 = Ka. 

This length scale a sets the natural length scale for the boundaries. We will assume 
that the boundary conditions are compatible, and further, the asymptotic shift is on 
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the scale aa^, i.e., 6 = Aaa"^, where A stays order 1 as we change the parameters of 
the problem. 

Since the linearized strain jyy is zero at the boundaries, we have = 1 — 
Outside \y\ < k, Wy = 0. Since the asymptotic shift is 6, we get w^(±L, y) —y = 
for \y\ > k. 

The elastic energy penalizes the square of the curvature. Consequently, all finite 
energy configurations of the sheet have a continuous tangent plane a.e. This yields 
the matching condition 

/ Vw^ \ _ f cos 2a —sin 2a \ / Vw^ \ 
\ Vv- J ~ \ sin2a cos2a J \ Vv^ J 

at y = 0, where V denotes the 2 dimensional gradient. Since (P) holds for all x at 
y = 0, it automatically implies that 

f '^x \ — ( '■^'^^ 2*-^ ~ ^^'^ '^^ \ ( \ 
\ '^x J ~ \ sin2a cos2a / V / 

We therefore have an independent matching condition 

(3) ('"^] = ( '"^o^"" )(''t) y = o 

^ ' \ Vy J \^ sm2a cos 2a J \ J 

For small a, we can linearize the matching conditions, and to first order in a, we 
obtain 

w~(x,0) = w^{x,Q) 
Wy{x,0) = Wy{x,0) — 2a 

(4) v-{x,0) = v+{x,0) + 2aw+{x,0) 

If we assume that Wy{x,0) ~ 0(a), then for small strains, Vy{x,0) ~ 1, and the 
matching condition for Vy is satisfied as an identity to first order in a. 

Since a is the natural length scale at the boundaries, Wy ~ 0(a) near y = 0,f ^ = 
?/± |5 and Vy +^{WyY — 1 = 0, we define the boundary conditions in terms of scaling 
functions - 

wf2is) = aaV2^f2 (-) ^ ^i2i-'^) = a^a^f^ (-) + s± aa^ — . 

' \a/ ' ' \a/ 2 

The scaling functions ($i'2('7)5 ^i'2(^)) normalized (have curvature 1) 
r ^\^tydri + f [$ r'Trfr/ = 1, for z = 1, 2, 

Jo J-L'/a 

are strain free, 

[$f]2 + ^f = 0, for 2 = 1,2, 
and satisfy the matching conditions 

$+(0) = $r(0) for 2 = 1,2, 
^r'(O) = $+'(0)-v^ for 2 = 1,2, 
^r(O) = ^+(0) + 2v^$+(0) + A fori = 1,2, 
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Finally, we also have the conditions v = ±|5, u = w = Oa.ty = ±L', where 6 is 
the asymptotic shift determined by the (compatible) frame boundaries at x = ±L. 
We will explore the relation between the frame and the value of A (and 

consequently the value 6) now. Using = for \y\ > k, and the definition of A, 
we see that 

^f{n) = ^f{ri) = for \r]\>K = k/a, i = 1,2. 
Using 'jyy = 0, we get 

^^fiv) = - [^f']'^dr] for ?7 > 0( resp. ?7 < 0). 

J±L'/a 

Using this relation in the matching condition, we obtain 



A 



2v^$+(0)+r mTdv+ r\'^t' 

J-L'/a Jo 



This gives the explicit relation between A and the out of plane displacement of a 
no-stretching profile. 

2.4. Rescalings. From Lobkovsky's results [2H] and the analysis in Re. [20]; we know 
that the dominant energy balance is between the longitudinal curvature w^y and the 
transverse strain w^. The natural length scale in X is L, and as we saw above, the 
natural scale for Wy \s a, the bending angle. 

Let Lq denote the natural scale for the quantity q. For example, we have = L. 
The above considerations lead to the conclusions 



Ly 



2 2 
(T W ~ 

yy 



r4 



LyY, 



This motivates the introduction of rescaled coordinates and displacements by 

x = LX, y = a'/'L'/'a~'/'Y 

and 



u 



y±- + ai/3L2/3a5/3y± 
^ 2 



X 



Here 6 is the asymptotic shift, and is determined by the frame boundaries. Since 
a, L, X, y, u, v, w, 6 all have dimensions of a length, and a is dimensionless, it is clear 
that the rescaled quantities X,Y,U,V^ are all dimensionless. Note that these 
rescalings are different from the rescalings in Ref . . 

With these rescalings, the dimensionless energy = a^^l'^L^^I^a^'^ I^X^ is given 

by 
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(5) 



+ e-^l\VY + Wlf\ + + 2^!^Wly + e^l^'Wlx] dXdY 



where we have to use the appropriate V"^ and W"^. We have suppressed the super- 
scripts ± on F and W for clarity, e = a /{La) is the natural small dimensionless 
parameter in the problem. We think of e as a dimensionless thickness. 

Our quest for rigorous scaling results for the energy I reduces to the following - 
Show that the rescaled energy 7 = 7+ + 7~, of a minimizer {U*,V*, W*), is bounded 
above and below by positive constants uniform in the dimensionless thickness param- 
eter e, as e — 0. 

With this rescaling, we have the matching conditions 

W+{X,0) = w-(x,o) 

W+{X,0) = Wy{X,0) + V2 
V-{X,0) = V+{X,0) + 2V2W+{X,0) + — 

Ly 

(6) = V+{X,0) + 2V2W+{X,0) + AA 

and the boundary conditions 

W^{-1, Y) = A^t (j) , U{-1, Y) = 0, V^{-1, Y) = (j 



U{1,Y) = 0, VHl,Y) = A^^ (^j 

(7) W^{X, ±L'/Ly) = 0, U{X, ±L'/Ly) = 0, V^{X, ±L' / Ly) = 0, 

where A — a/Ly is the nondimensionalized inverse curvature, and ^f,^f are the 
scaling functions from above. 

3. Existence of a minimizer 

Prom this point forward, we will use c, C, C , Ci, C2, etc. to denote constants whose 
precise numerical values are not important. These constants can change from one 
line to the next. By doing this, we can focus on the scalings of the various quantities, 
without worry about the numerical values of the constants in the scaling relations. If 
these constants depend on a parameter we show this dependence by writing C{q) 
or Cq. Also, we will suppress the superscripts ± whenever this does not cause any 
confusion. 

Lemma 3.1. u, : = [— L, L] x [0, ±L'] R are 77^ functions such that u — Q 
at X = ±L and m = at y = ±L'. The functional 

1/_L nO // 2|/-\2 



J{u,v+,v-) = JJ^^K + iv^)' + -{v: + Uy)' + JJ^ K + {Vy)' + -{v; +Uy)' 
is coercive, in the sense 3c, C > such that 
J{u,v^,v-) > c{\\Du\\l2 + \\Dv+\\l2 + \\Dv-\\l2)-C / [v+{x,0) - v-{x,0)]'^dx 



Proof. We first consider the integral over S~^. Using the boundary conditions, and 
repeated integration by parts yields, 

Uy{x,y)v^{x,y)dxdy = — J u{x,0)v^{x,0)dx 

uix,y)v^y{x,y)dxdy 
u{x, 0)v^{x, 0)dx 

-L 

(8) + jj u^{x,y)v^{x,y)dxdy 

Adding a similar result for S~ , we obtain 

1 /"^ 

J = -{\\Du\\l2 + \\Dv+\\l2 + \\Dv-\\l2)+ u{x,0)[v-{x,0) -v^{x,0)]dx 

^ J-L 



1 rr , , . ..o. . 1 
'4 



+^ j j ^{ux{x,y) + v:^{x,y)ydxdy + ^ jj {u^{x,y) + {x,y)fdxdy 



> ^iWDuWl^ + WDvm + WDv-Wl) 
(9) - / \eu'^{x,0) + Ce\v-{x,0) -v^{x,0)\ 



dx 



where ^ > can be arbitrarily small. By the trace theorem and the boundary 
conditions u = at a; = ±L and y = ±L', we have 

j u'^{x,0)dx < C\\Du\\l2. 

Using this inequality in the estimate (jHl) for a sufficiently small 6 we get 



J{u,v^,v-) > c{\\Du\\l2 + \\Dv+\\l2 + \\Dv-\\l2)-Ca^ 



[v^{x, 0) — {x, 0)]^dx 



and we can take c = 1/4 if we so choose. □ 

We are now in a position to prove the existence of a minimizer for the function X 
with 0" > 0, for our "no stretch" boundary conditions. 

Theorem 3.2. If a > 0, every minimizing sequence {uj,v^,w^) for the energy func- 
tional I, that satisfies the boundary and the matching conditions, has a subsequence 
that converges in H^{S) x H^{S) x H'^{S) to a global minimizer of 2. 

Proof. This conclusion follows easily from the direct method in the Calculus of vari- 
ations nil Emm. 

{uj, v^, wf) is a minimizing sequence for X that satisfies all the boundary conditions 
and the matching conditions at y = 0. For all j, we have the boundary conditions 
wf = v^a;a<l>^2(?//o) at x = ±L and wf = at y = ±L'. These boundary conditions, 
along with X > cr^||-D^w||^2, imply that Wj is a bounded sequence in H'^. 
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Consequently, up to extraction of a subsequence, wf w^* in H'^. The compact- 
ness of the embedding H'^ — > W^'"^ imphes that (w^)^, {WyY converge strongly to 
(w^*)^ and {w^*y respectively in L^. 

The standard trace theorems [3H] imply that 

w^*{±L, y) = y2aa$^;2 {£) , w^*{x, ±L') = 0. 

Consequently, {uj,vf,w^*) satisfy the no-stretch boundary conditions for all j. In 
general, they do not satisfy the matching conditions at y = 0. Rather, we have the 
relation 

v+{x, 0) - vr{x, 0) = 2awf{x, 0). 

We set We will henceforth drop the subscript j where it wont cause any 

confusion. Since {uj,vf ^w^*) satisfy the no-stretch boundary conditions, -S = at 
X = ±L and at y = ±L'. Also, from the matching conditions for f^, we have 



v^{x, 0) — (x, 0) = 2aw^{x, 0). 



Since llD^ii?^!! is bounded, and 



/ (w+(±L, CaV, / {w^{x,L')fdx = Q, 

Jo J-L 

it follows from the trace theorem that w^{x, 0) is bounded in L^. 

With this definition of it, and suppressing the subscripts j and the superscripts ±, 
we have 



I{u,v,w) > j (^u^ + ^^lj + 



1 \ ^ 1 _ 1 

' + Uy + W^Wyf + {Vy + ^W^^'^ 



dxdy 



Given any e > 0, up to extraction of a further subsequence, we get 

1> jj {u.f + {vyf + \{v. + uyf - C{\\w*\\l,,, + e)(||D^|| + \\Dv\\). 

The argument from above shows that f^!" — = 2aw^ is bounded in L^. Combining 
this with the previous lemma, we see that \\Du\\ + ||-Df || is bounded in L^. 

u satisfies the boundary conditions tt = at x = ±L, y = ±L'. v satisfies the 
boundary conditions 

v^{±L, y) = y±^- + a'a^f^, (^) v^{x, ±L') = ±L' ± |. 

By the boundedness of ||-Df ||, it follows that 0) exist in the sense of traces, and 
further are in L?[dx). Consequently both u and are bounded in L'^[dS^). Combin- 
ing this with the boundedness of \\Du\\ + ||-Df ||, it follows that a further subsequence 
{uj,Vj) converges weakly to {u*,v*) in x H^, and consequently {uj,Vj) {u*,v*) 
where 

U* = U* + X 

Since Hq is weakly closed, it follows that u* and f ^* satisfy the boundary conditions 
in the sense of traces (See jTHj for a detailed argument). We now observe that the 
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functional X is convex in Du, Dv and D'^w, and consequently weakly lower semi- 
continuous on X X H^. This implies that (u*, f*, w*) is a minimizer for I. □ 

The existence of a minimizer opens the door to a direct analysis of the Euler- 
Lagrange equations for the energy functional X. We will pursue this approach else- 
where. In this paper, we restrict ourselves to arguments that do not involve forming 
the first variation of X. 



4. Lower Bound 



In this section, we prove a lower bound for the linearized Elastic energy X in Eq. 0, 
by proving a corresponding result in terms of the scaled energy / in Eq. (0). In the 
remainder of this section, we will mostly work with the half sheet , although 
the same arguments also hold on S~ . With this understanding, we will drop the 
superscripts ±. 

As we show in jS^, it follows from the boundary condition U = OatX = ±l, that 



(10) 



> 







1 N 2 

2 



-1 



W^ydX 



-1 



dY 



So, it suffices to prove a lower bound for the functional 



E{W) 







WldX 



-1 



W^ydX 



1 



dY 



As in 
(11) 



we let E}, and Eg denote the quantities 

fOO 1 



Eu 



W^ydXdY, 



E, 







1 N 2 

WldX ) dY. 



which we will call the (lower bounds for) the bending and stretching energies respec- 
tively. 

For every X, we define 

-1 



P{X) 



W^YiX,Y)dY 



p{X) is a "local" (in X) measure of the bending energy, and [p(X)]~^ can be thought 
of as the bending energy density in X that is obtained by integrating out the Y 
dependence. 

What we will see below it that p(X) is the natural length scale associated with the 
ridge as a function of X, viz., Wo{X) ~ p(X), and the bending energy density in Y 
decays rapidly for Y/p{X) ^ 1 (See Fig. H]). 

Before we begin the proof of the lower bound, we prove the following elementary, 
but very useful result. 

12 



Lemma 4.1. f e H\ J {f'f = p-^ < oo, /(O) = /o, /'(O) = (3. Then we have 



/ f'^(ri)d'n > max (1 
Jo ^ ' ' - z<y, ee(o,i) ^ 



Proof. By the Sobolev Embedding theorem, / is a function and 

f{0 = h + K+ t f"ivm-v)dv- 

Jo 

Defining T(0 by T{^) = f"{r]){^ - r])dri, we have 

Jo Jo '~>P 

and integrating this equation in ^ yields 

If t{Z) = p{ri)dri, t is nondecreasing in Z. Using f{Z) = fo + (3Z + T{Z) and 
the elementary inequality 

|a + 6p > (l-^)|a|2-l^|6|2, 



9 



for all < < 1, we see that 



riZ) > {1-6) ifo + f3vrdv 
Jo 



2.. 1-^ f\2 



> (I 



> a 



f^Z + f3foZ^ + ^Z']- 



3 



T {ri)dri 



1-6 



6 



12p 



The result follows by observing that r is nondecreasing, and optimizing the choice of 
6. □ 

Theorem 4.2. I{u,v,w) is as defined in Eq. For all u G H^,v^ G and 
G n W^'"^ satisfying the no-stretch boundary conditions 

and the matching condition Wy{x, 0) = Wy{x, 0) + 2a, we have the lower bound 



I{u, V ,w ) 



> min (^C^a'/'a'/'L'/^C2a'a'^^ . 



Remark. Note that we do not need all the matching conditions in Eq. (^. 

13 



Remark. The form of the lower bound gives a crossover scale La 



~ (J 



L/3^2/3^-1/3 



Ly for a. This scale is the same as the one we obtained in our earlier rescalings. 

We will prove the theorem by doing the cases a -C Lq and a ^ La separately. The 
result for a <C is obtained by proving the scaled version of the statement, viz., 

I >Ci ioT A< A*, 

for a constant A* that will be determined below. The case A > A* is much easier, 
and follows immediately as a corollary. 

Remark. In earlier work jSE], we proved the same scaling result for a <C La, but 
with extra assumptions on the behavior of v and w at y = 0. We do not know a 
priori that these assumptions are satisfied for a real crumpled sheet. In this theorem, 
we have removed these hypothesis, and this result is directly applicable to crumpled 
sheets. 



We now begin our proof of the theorem for a La- As in j2SI, the idea behind 
the proof is to show that the stretching energy can be bounded from below by a 
negative power of the bending energy so that the total energy Eg + Ei, tends to 
+00 as — > and Ei, — >• 00. This ensures the existence of a positive lower bound 
for E (and consequently also for /). 

We set (3^{X) = Wp{X, 0). The matching condition therefore is = + 

Before we prove the theorem, we collect a few useful results in the following 
lemmas. 



Lemma 4.3. 



Y 



1 Jo 



Wi dYdX > C max 

Z<Y 



l3{XfdX - 2Z^Eb 



C'A^ 



Proof. Using Lemma l4.ll with /(^) = W{X,^), taking 6 = 1/2 and integrating the 
result in X, we see that 



1 rY 



1 Jo 



W'^ dVdX > — max 
~ 24 z<Y 



j (3{X fdX - 2Z^Eb 



The Poincare inequality now yields 

•1 i-Y rl pY 



/I /.y pi pY pY 

/ WldYdX>C / W^dYdX-C [W^{-l,i) + W^{l,i)\di. 
-1 Jo J~l JO JO 



The result follows from the observation 



pY poo 

/ [W\±l, 0] di<A^ ^l2{v)dri < C A^ 
JO Jo 



□ 



Our proof is based on demonstrating that a small bending energy E^ will lead to 
a large stretching energy. For A < fi, this idea is quantified by the following lemma. 
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1 /2AEfe^ ^ 



Lemma 4.4. Let B = j^^j3{XfdX, and let 

^ B\ B 

There is a constant /i* > such that, if < yU*, the stretching energy i^^ and the 
total energy satisfy lower bounds 

CB'' 



El 



and 



E>Eo = i5CB 



7N1/6 



Proof. By Jensen's inequality, we have 
E.= I U I 



WldX 1 dY > 



1 

2 Jo 



1 dY > 



1 

2F 



Y rl -I 2 



2 vj_i 
Lemma (4.31 now implies that 

> max 2-1/2 ^c{BY'/^ - 2E,Y'''^) - C'A'Y-^I^] . 

Setting i?y^/2 = 2y^/2ii^ft, we deduce that a characteristic scale Y for F is given by 

B 



Y 



2E. 



Rescaling Y in terms of Y, we obtain 



C57/2 



5/2 



I) 



5/2 



where fj, is as defined above, i.e. 



1 /2AEA' 



Observe that has a positive maximum a.t z = 5/7 > 0. The lower bound for 

the stretching energy Eg follows by continuity of the function — C'fiz ' 

with respect to /x at ^ = 5/7. 
Minimizing Eg + E^, we see that 

E = Es + Eb>Eo = {5CBy/'' 

□ 

We can now prove the theorem. 

Proof. Set A* = {n*f/^/{2Eo), where /i* and E^ are as in Lemma [4.41 
The precise statement we will prove is 

E > Eomm{l,A*/A) 

If Ej, > Eq, there is nothing to prove. Therefore, we can assume that E^ < Eq. It 
follows that Eb < Eq for each half-sheet. 
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We first consider the case A < A*. As we argued before, tlie slieet lias a well defined 
tangent vector at ?/ = 0, and this gives the matching condition (3+{X) = P-{X) + \/2. 
This, along with the convexity of the map (3{X) ^ B = J\f3^{X)dX, implies that 
B+ + B- is minimized when P+{X) = -p'{X) = ^. Therefore, W.L.O.G. B+ > 1. 
For the half sheet the half-sheet y > 0, we have - 

A<A*,E<Eo and B>2 imphes that 

1 f2EbA^ ^ 



Lemma (4.41 now implies that E > Eq, for the half-sheet y >0, and consequently for 
the whole sheet. 

If y4 > A*, we still obtain the conclusion E > Eq hj the preceding argument if Eh 
is so small that 

1 /'2EhA^ ^ 



Therefore, we only need to consider the case fi > fi*. W.L.O.G B{= B^) > 1, so 
that fi > fi* implies 

and this gives the desired conclusion. □ 

5. Upper bounds 

Our goal is to obtain upper bounds for the functional / that scale in the same way 
as the lower bound from the previous section as a function of the nondimensional 
parameters in the problem, viz. e, a and A. This will show that we have captured 
the optimal scaling behavior of the elastic energy for a single ridge in a crumpled 
sheet. 

In particular, we want an upper bound that is a constant (independent of e, a and 
A) a A < A*, and an upper bound that scales as 1/A for A > A*. Also, we want 
upper bounds that are independent of e and a. 

The existence of such upper bounds can be motivated as follows. The energy in 
the rescaled variables is given by 

r r r e-^/s , i 

I{U,V,W) = J J {Ux + W^y + ^{Vx + UY + 2WxWYy + e-'/^VY + W^f 

+ [W^Y + 2e^/^W^Y + e^^^W^lx] dXdY. 

We would like to show the existence of (f/, V, W) satisfying the boundary conditions 
such that I{U, V, W) < C{A) < oo uniformly in e for A > 0. The idea behind the 
construction of an appropriate {U, V, W) is as follows. We first pick a smooth W 
satisfying all the boundary conditions. For this W, we will pick an V such that 
Vy = —Wy- This equation can (we hope) be solved for every X, along with the 
appropriate boundary conditions V{X,Y) — > as F — ±oo. Once we have V, we 

16 



determine U hy Uy = —Vx — 2WxWy, again with the appropriate initial condition 
for U. With such a choice for U, V and W, the energy becomes 



j I [{Ux + Wlf + W^y + + e'l^Wlx] dXdY, 

and since U, V and W are assumed smooth, it easily follows that there is a finite 
upper bound, uniform in e as e — *• 0. Of course, we are not guaranteed that we have 
the right dependence on A. Also, we are not guaranteed to get the right asymptotic 
shifts in V^{X,.). 

In the remainder of this section, we will deduce the upper bound by using ideas 
similar to the simple argument from above to explicitly construct smooth functions 
(m, f w"^) satisfying all the boundary conditions. With these functions, we can show 



X(m, V, w) < min 



Ca'/^L'^W/' + C'aW log+ (^) , C'a^a 



a 



where log^ x = max(logx, 0). This is not exactly the scaling that we obtained for the 
lower bounds. In particular, the upper bound indicates that we are missing some of 
the relevant physics in our lower bound if a <^ crexp(— e"^/'^), i.e A <^ e^/^ exp(e~^/'^). 



5.1. Self similar test solutions. In this section, we show that, for identical, no- 
stretch boundary conditions at x = ±L, with zero asymptotic shift (A = 0), we can 
construct "self-similar" test solutions that yield the "correct" upper bound. 

Let $^ and be as in the definition of the boundary conditions, so that $^(77), \E'^(?7) 
are smooth functions, that are supported 'm.\ri\ < K. 

We will choose W and V in the following "self-similar" form 

wHx. V) ^ Pix)^- (^) , vHx. V) ^ Pix)^- (^) , 

where p{X) is a function smooth that will be chosen later satisfying p{X) > for all 
X. The boundary conditions at X = ±1 require that p(l) = p(— 1) = A. 

Let T] denote Y/p{X), so that W^{X,Y) = p{X)^^{r]). From this we obtain. 

Wp{X, Y) = $±'(r/), V^{X, Y) = vl/±'(r^) 



From the matching conditions on $ and ^ at 77 = 0, it is clear that the functions 
and from above satisfy the appropriate matching conditions © at y = 0, 
only if A = 0. 

Differentiating in X we obtain. 
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Differentiating once more, we get 



W$y{X,Y) = 
W^y{X,Y) = 



——^^"in) 
p{xf 



From tlie no-stretcli boundary condition + ^ = 0, it follows that Vy + 



0. 



U^{X,Y)=p'{X)p{X)E^^ ^ 



We would also like Up = —V^ — 2WxWp . Prom the above scalings, we see that 
+ 2WpWp = p\X) ^±(77) - r]^^'{r]) + 2^^' [r]) (^^{r]) - r)^^\r))^ 
In order that Up have this scaling behavior, we will set 

If we choose S such that 

S±' = - [^±(77) - v'^^'iv) + 2^^'(?7) {^^{v) - V^^'iv)) 

Then, we will have Up = -V^ - 2WpWp. 

In this approach, we have first order ODEs for S^, with boundary conditions 
5^(77) ^ as 1] ^ ±00, and these ODEs can be solved (in principle) to yield the 
functions S^. These functions are also required to satisfy the matching condition 
U+{X,0) = U-{X,0) for all X, i.e, the condition S+(0) = S-(0). 

Since the unique solutions for the ODEs determining 2=*= are also supported in 
\v\ < K. 

Using the no-stretch condition + ^ = in the ODEs for 2=*=, we get 

= - Uf^if}) + r/^^'(r/) + 2$^' (77) $^(77) 

Integrating the above equation, using the fact that $^ and \E'^ are supported in 
I77I < we deduce that is supported in [—K,K] and S='=(0) = [$=''(0)]^. The 
matching condition for now implies that S+(0) — H~(0). Consequently, 



U^{X,Y)^p\X)p{X)'E 



\p{x)) 



satisfies the boundary conditions U^ — > as y — > ±00, and the matching conditions 
U+{X, 0) = U-{X, 0) = 0. Wc also require that [/ = at X = ±1. since p(±l) = A, 
we will now require that p'(±l) = 0. 

We will henceforth restrict ourselves to considering the half-sheet F > 0, since the 
same arguments will also apply to the half-sheet ^ < 0, and we can now drop the 
subscripts ±. The above procedure yields an appropriate test configuration ([/, V, W) 
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for boundary conditions that are identical at x = ±L and satisfy A = 0. We will 
henceforth refer to this situation as the self-similar case. In the remainder of this 
section, we will consider this special case, and we will consider the general case in 
Sec.lO 

From the above arguments we see that, for appropriate boundary conditions, it is 
indeed possible to choose {U, V, W) in the self-similar form 

W{X,Y) = p(X)$ ' ^ ^ 
V{X,Y) = p(X)vl/ 



PiX)J 



f/(x,r) = p'{x)p{x)E 



„ f Y \ 



such that Vy + Wy = 0, Uy + Vx + 2WxWy = 0. With these choices, we have 

I{U,V.W) = J J {{Ux^Wlf ^WlY^2e^l^Wly^t^l^Wlx\ dXdY. 

A straightforward calculation allows us to estimate the various terms in this expres- 
sion. We obtain. 



{Ux + WlfdXdY < C j {[p"{X)p{X)f + [p'{X)f) p{X)dX. 

1 



W^ydXdY < C I -^dX 



WlydXdY < C I i^^t^rfx 



(12) j I WlxdXdY < cj (^[p"(X)]V(X) + ^^||y^) dX 
where C is a constant that only depends on $, and 

5.2. Construction of the upper bound. In order to prove the claimed upper 
bound, we need to show the existence of a smooth p{X) such that p(±l) = A, 
p'(±l) = 0, and all the terms in / are bounded uniformly in e as e — * with the 
appropriate dependence A. 

Lobkovsky's [2H1 analysis motivates the choice p(X) ~ (1 — X)^/^ near X = 1. For 
this choice however, the contribution of is given by 

[pWpW-i^^-(i-^n 

and is not integrable near X = 1. Therefore the choice p(X) ~ (1 - |X|)2/3 will not 
yield an upper bound for /. 

In our analysis of the lower bound, we ignored the contribution of Wxx to the 
energy, and obtained results that agree with Lobkovsky's boundary layer analysis. 
This suggests that the scaling p(X) ~ (1 — |X|)^/^ might still be appropriate in 
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regions where the contribution of the Wxx and the Wxy terms are small. However, 
close to the boundaries near X = ±1, the dominant energy balance is different, and 
we need to modify the behavior of p to account for this. 

For small X, we expect that e'^^^Wxx ^ Wyy is the leading order balance for the 
energy m /. Since Y ~ p(X), it follows that p{X) ~ e-V3(i _ x) near X — 1, 
and similarly p{X) ~ e~^/^(l + X) near X = —1. Note that these behaviors match 
p{X) ~ {l-Xy/^ (respectively (l+X)^/^)) when 1-X ~ e (respectively 1+X ~ e). 

We also have the boundary conditions p(±l) = A and p'(±l) = 0. This suggests 
p{X) K, A ioi {I - \X\) < Ae^l^ if A -C e^/^ Therefore we wiU choose p{X) with the 
following behavior. 

• In the case A < e^/^ < 1, 



-V3(l_ |X|) 



(1 - X2)2/3 

In the case e^/^ < A < 1 



(l-|^l)<^eV3 

AeV3 ^ (1 _ |x|) < e 
e< (1- |X|) ~ 1 



(1-I^l)<^ 



[1-X'^fl^ e<(l-|X| 

• In the case A ~ 1, we set piX) — A. 

We will now make the above considerations precise. Let Lp G 
nonnegative, non-increasing function, that is identically one on 
on [2, oo). An example of such a function is illustrated in Fig. Also, (p will denote 
the complementary function 1 — y?. 

Lemma 5.1. For A > and e > 0, let 



be a smooth, 
oo, 1/2] and zero 



9{z) 
h{z) 



where x denotes the characteristic function. Then is a smooth function. Also, there 
exist constants c and C such that g satisfies the following inequalities \/z 



9{z) 
9{z) 



> 
< 

< 



c [^X[o,2A6V3] + h{z)] 

C [^X[0,2Aei/3] + K^)] 



Ch{z) 
g'\z) < Ch{z) 



1 

1 + - 

z 



1 + - 
z^ 



Proof. We begin with an elementary observation. Near z ~ Ae^/^, the functions 
g = A and g = ze~^l^ are comparable, ie, there exist constants 6', independent of 
e. A, such that OA < ze~'^/^ < OA for ^Ae^/^ < z < 2Ae^/^. Similarly, near 2; ~ e, we 
have ^ze-V3 < ^2/3 < Qze'^/s for le < 2e. 
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This observation implies that g{z) > c [^X[o,2Aei/3] + h{z)'j . The inequahty g{z) < 
C [^X[o,2Aei/3] + h{z)\ is elementary and follows from the boundedness of 
We also observe that, for alH > and n = 1, 2, 3, . . ., 



dz^'^yi 



Cn 



and the same inequality also holds for the complementary function ip. We also have 
the elementary inequality 



dz'^ 



< — ■ , z > 



Differentiating g twice, using all of the above observations, and recognizing that 
z~^ < 1 + we get the inequalities 



\g\z)\ < Ch{z) 
\g"{z)\ < Ch{z) 



z 



1 + - 

z^ 



This proves the lemma. 



□ 



Lemma 5.2. For boundary conditions that support a self-similar test function, we 
have an upper bound 



/([/*, < C 



1 + e^/^ log"* 



,2/3 



Proof. We set p{X) = g{l — X^), where g is as defined in lemma 15.11 Note that 
^(0) = A, c/'(0) = implies p(±l) = A,p'(±l) = 0. 
We begin with a few observations. 

(1) Our definition makes p an even function of X, and by lemma IKTTI p is smooth. 

(2) z'{X) = -2X and z"{X) = -2 are bounded for X e [-1,1], so that the 
inequalities in lemma 15.11 also hold in terms of p with the understanding 
z=l-X\ 

(3) If A < e^/^ < 1, the sets [0, Ae^^^), [Ae^^^, e) and [e, 1] give a disjoint partition 
of [0, 1]. If A > e2/3, the sets [0, Ae^/^), and [e, 1] cover [0, 1]. 

We can now estimate the various terms in ()12|) using the results of lemma 15. H and 
the observations from above. 
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We first estimate the {Ux+W^y term. Using {l+z'^ < Cil+z'^ < C{l+z-^), 
we obtain 

J {[p"{x)p{x)r+[p'{x)r)p{x)dx < cj\h{z)r(i+^yz, 

° i: 



z^ I 1 



dz 



dz^ 



The Wyy term yields 



dz 



+ C 



J^Ae^lz Ze V3 



dz 



< C 



3g 



,2/3 



+ 1+1 



We will now estimate the Wxx term, since a bound for the Wxy term can be obtained 
from the bounds for the Wyy and the Wxx terms. We have 

[p"{x)fp{x)+''p^yx < c j\h{z)r (^1 + ^yz 



< Ce-^ f z^ (l + X]dz 



+c 



Using 



[p'{x)r [p'{x)Y- 



p\x) p{x) p{xy 

and the Cauchy-Schwarz inequality, the above estimates yield 



/ 



P{X) 



■dX < c 



,2/3 



A 



+ 1 + e- 



-2/3 



We have thus bounded the Wxv term. Using these estimates in 

I{U, V.W) ^ JJ [{Ux + W^y + W^Y + '^e^'^W^Y + e^^^W^lx] dXdY, 
we see that, for A < e^/^ the self-similar test function has 



KC 



1 + log 
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,2/3 



If 6^/3 < A, the sets [0, Ae^/^) and [e, 1] cover [0, 1]. In this case, the same analysis as 
above gives I < C. Combining this with the above result proves the lemma. □ 



Setting p{X) = A, a. direct calculation using (fT^ . shows that / < C/A. Combining 
this result with the lemma and "unsealing" these results, we get 



Theorem 5.3. // the no stretch boundary conditions are identical at x 
have A = 0, we have the upper bound 



±L, and 



I{U\V\W*) < Cmin 



1 + e^'^ log^ 



,2/3 

~A 



A 



5.3. Non-self similar test solutions. We will now consider upper bounds for the 
general case, i. e, for the situation where the boundary conditions at a; = ±L are not 
identical and/or A 7^ 0. 

The strategy of the proof will be the following: 

(1) We introduce boundary layers (in x) near the boundaries x = ±L of width 
X = b. In these boundary layers, we connect the boundary condition at 
X = ±L with identical profiles that have A = 0, at a distance a from the 
boundaries. 

(2) In the region, L — |x| > b, \y\ < KLy, we use our self similar construction 
from the last section. 

(3) In the remaining region, we introduce a small, uniform strain ■jyy ~ S/L', to 
get the appropriate asymptotic shift S. 

The various regions are illustrated in Fig. |21 The idea behind this construction is 
from matched asymptotic expansions. The self-similar solutions from the last section 
play the role of the outer solutions in x, but the inner solutions in y\ 

Before we rigorously construct a test solution that gives the upper bound, we 
present a heuristic scaling argument that motivates our choices for the length scales 
b and / in our test solution. 

The boundary conditions at a; = ±L are given by no-stretch profiles, whose (cur- 
vature) length scale is a. We assume that we can find a no-stretch profile with zero 
asymptotic shift, and curvature a~^. This will be the profile of the test solution at 
X = ±(L — b). Also, we set u = x in these layers, so that Ux = l,Uy = 0. 

In the boundary layer(s) L — b < \x\ < L, w{x,y) is supported in \y\ < Ka. Also, 
w ~ 0{aa) so that 

aa aa a 

b V a 

and all of these derivatives are supported in \y\ < Ka. This gives '-yxx ~ "U^^ ~ o?c? jl? 

If we assume that there is no stretching in the y direction, so that ^yy = 0, we see 
that f y ~ 1 — ^Wy. Consequently, we see that the asymptotic shift 5 is a function of 
X, and by integrating Wy we see that S{x) ~ 0{a^a), and this gives ~ ~ a^a/b. 
Also, since this Vx arises from the difference in the asymptotic shifts, it is supported 
in \y\ < I, and not only in \y\ < a, as for Wy. This gives 'jxy ^ Vx ^ a^a/b. 
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Figure 3. A schematic illustration of the various regions, and their sizes. 

Putting all of this together, the elastic energy in the two boundary layers compris- 
ing Region I is 

~ llx(A + llylh + CT^ [wl^ + wly) ah 

a^a^ a'^a?l a'^a'^b a'^a^a? 
+ —— + + 



h a 63 

In Region II, we use the self similar construction from above, i.e., the dominant 
energy balance is between the curvature Wyy and the strain '^^x- Following the scaling 
argument in Ref. |^, we get 



£ 



II 



L3 / 



In region III, we are connecting v = y (zero asymptotic shift) with v = y ± 6. We 
can do this with a profile that has a uniform strain 'jyy ~ S/{L' — I). Consequently, 



^L'-IJ ' ' L'-l 

Note that Sjjj is completely determined by the boundary conditions, and is indepen- 
dent of any choice we make for b and / provided that I ^ V . 

For any given a, a, L, L', £ = £i + £ii + £in diverges as 6, / ^ 0, oo independently. 
Consequently there are optimal finite, nonzero choice b = b* and 1 = 1*. 
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There are two scaling regimes of interest. In physically realistic situations, cr < a. 
li a — Ca, with C staying as cr —> 0, optimizing I and b gives 

I* ~ a-^/^a'/'L''/^ b* ~ a'/^e^/'L^/^ 

The energies in the three regions are 

For a <ti L, L', the energy is asymptotically negligible. £j ~ a^^'^Sjj, so the two 
energies scale in the same way for a ~ 0(1), but the energy in the boundary layers 
is asymptotically negligible in the small angle limit. 

We can also consider the situation a <^ a. In this case, the natural scaling regime 
is a ~ Ca^l^L~'^l'^ ^ with C staying as cr — > 0. Optimizing I and b gives 

r ~ Q;-V3^V3^2/3^ ^* ^ ,5/3^-2/3 ^ ^ 

The energies in the various regions are 

El ~ o?a\ En ~ a'/V^/^^^/^^ ~ aV^ (^) 

For a <^L^ L', the energies Ei and £^7// are asymptotically negligible, and the energy 
is determined, essentially by the the self-similar solution in Region II, in the limit 
e = (a/La) ^ 0. 

We will now use these scaling results as motivation, and rigorously construct test 
solutions that give the appropriate upper bound in situations where the boundary 
conditions are non-identical, or have nonzero asymptotic shift. Before we begin our 
construction, we first show that there exist no-stretch profiles with zero asymptotic 
shift. 

Lemma 5.4. 3^0, such that > Kq^ there exist smooth 0, il) supported on [0, K\ 
such that 

-K 



and 



Proof. Let ^ be a smooth nonnegative function such that < ^ < 1, ^(3:^) = 1 for 
X < —1 and ({x) =0 for x > 0. Let ci7 be a smooth function supported in [—1,0], 
that is odd about — |, and is not identically zero. 
We set q{x) — ({x) + qw{x). The map 

g^A(g) = y"\(o [i-^mdi 

is clearly continuous. The reason for calling this map A will become clear below. 
We have ^ 

A(o) = j [i-c(e)]t^e>o 
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since < C < 1- Also, 



< 



aod^ 



-1 



-1 



where we have used |C| < 1 and w is not identically zero in passing to the last line. 
Since A(0) > 0, and A(g) -oo as g — > -oo, 3g* < oo such that A(g*) = 0. 
We set ^ = ^ + q*w, and define / by 



-1 



Let = 21. For a given K, we will set 



K 



Then, (p'{x) = q {^), so that (f)'{x) = 1 for x < I < K - 1. Also, 



[(p" (x)]'^ dx 



1 



K r 



x-K 



dx 



1 
2' 



Since is supported on [0, K], it follows that 

-K 



(t)'{x)dx. 



Since = 1 for < x < — /, it follows that 



[<p\v)fdv+m 



K 



K-l 



[{<P\n)f - 0'(^)] dr, 



We now set 



ip{x) 



K 



A{q*)l = 0. 



□ 



We can extend the functions 0, to a no-stretch profile on M by setting 

^+[r]) = ^'{-1]) = (f){r]), ^+{r]) = -^-{-7]) =ilj{r]), r]e[0,K] 

and $^ = = otherwise. This profile satisfies the matching conditions and has 
zero asymptotic shift. 

We now have to prove that, by introducing a thin boundary layer near the bound- 
aries, we can "connect" the prescribed boundary conditions to the profile we con- 
structed in lemma 15. 4[ without incurring a large energy penalty. We begin with the 
following lemma which estimates norms of the first derivatives of a no stretch profile 
($^, \Ef^) in terms of the curvature of the profile, and it's support. 
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Lemma 5.5. ($'^, '^^) is a normalized no-stretch profile with support K, i.e., \E'^' + 
[$±']2 = 0, $±(r^) = ^±(ry) = if |r/| > K, 



$-'(0) = $+'(0) - V2 



and 



(?7) (i?7 
It then follows K >l/2 and 



<l>""(r7) (ir] = 1. 



(13) 



1 



K 



< C 



where C is a universal constant (independent of K). A similar result also hold for 

Proof. This follows immediately from the Sobolev embedding theorem in M} \n\ l38] . 
Since $^ {K) = 0, we have 



|$+'(r/)| 



< 



1/2 



< ^/K-r] 

Using $"'(0) = $+'(0)- v^, it follows that 2/^ > ^2, which implies that K > 1/2. 

The remaining inequalities follows from integrating |$+'(?7)| < -/K~^^ and using 
$+(ir) = 0, and is: > 1/2 so that > ij^3/2_ □ 

Lemma 5.6. $^''2 and '^^2 smooth on and supported on [0, ±K]. Further, 
for i = 1, 2, + [$f ]2 = 0, and 



$+"(77) dr]+ $ "(77) dr] = l. 



Also, $^''2 and \E']''2; satisfy the matching conditions (jHI), where the asymptotic shifts 
of the two profiles are given by 



1 2 



$-'(//) cir/ + v^<l>,^(0)+ / $+'(7/) rfr7 + y2$+(0) 



A,; 
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is a nonnegative, nondecreasing C°° function such that ({x) = for a; < and 
C(x) = 1 fo: 



> 1. Let 




UiX,Y) 




W^{X,Y) 








V^iX,Y) 


J±K 



c 



where 



A(X) = - 



/ [Wp{X,Y)fdY + Y±\/\{X) 

J±K ^ 



/•oo 

[WYiX,r])Ydr]+ / [W^{X,r])Ydr] + 2V2W'^iX,0) 
Jo 



Then, {U, V^, W^) satisfies the matching conditions at F = 0, and we have 



/ / {Ux + WlfdXdY < ^ 



M pB 



{Vx + UY + ^WxWyfdXdY < 



CM 



y=o J x=Q 



B 



pM pB 

/ / {Vy + W^f dXdY = 
Jy=o Jx=o 



M pB 

/ W^ydXdY < CB 

Y=0 Jx=0 
M pB ^ 

/ W'xxdXdY < — 
where C is a constant that only depends on K. 

Proof. Most of the inequahties follow by direct integration using the definition of 
{U, V^, W^), and using Lemma (5 .51 The only result that needs proof is 

r-M pB 



PIVI pH 

/ / {Vx + Uy + 2WxWYf < 



CM 



Since W{X, .) is supported in [0,K], and V^{X,Y) = Y± |A(X) for \Y\ > K, we 
have Vx = Ax- From the definition of A{X), we see that 

r-O 



d 
dX 



A{X) 



B 



Wy{X,v)[<l>,{7^)-<^^{v)]dV 



-K 



+ 



K 



W^{X, r7)[$+(r7) - $+(r/)]rfr/ + y2($2(0) - $1(0)) 



Estimating the integrals by the Cauchy-Schwarz inequality, and using Lemma 15.51 
we see that 



AvI < 



C_ 
B 



2-^ .C 



[^xiOYd^< 



B' 



28 



We also have, 

-^{yx + UY + 2WxWY) = Vxy + 2WxyWy + 2WxWyy 

= -^{Vy + W^) + 2WxWyy = 2WxWyy 
Using U = 0, and Vx = Ax, Wx = Wy = for Y > K, we get 

/ {Vx + UY + 2WxWYydr] = {M-K)A\+ / {Vx + Uy + 2WxWYfdT] 
Jo Jo 

pK 

< CMA\ + CK"^ I WlW^ydV- 
Jo 

using the Poincare inequahty. By the convexity of the map, 

-K 

WyYdri, 



W{X,.)^ [ 
Jo 



and the normahzation of ^^2^ follows that WyYdrj < 1. Lemma f5 . 51 along with 



implies that 

snp\Wx{X,7^)\ < 1. 

X,Y -D 

Combining this with the earlier estimate, we see that 



PIVI 

/ (Vx + Uy + 2WxWYfd7] < CMA\ + 
Jo 
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Integrating this inequality in X, and using M > K, we obtain the desired result. □ 

This Lemma provides a rigorous basis for our heuristic calculation for the energies 
of the three regions at the beginning of this section. Using this lemma, and the ideas 
form the heuristic calculation, we have 

Theorem 5.7. The minimizer of the elastic energy functional X in (0j subject to 
compatible, no-stretch boundary conditions 



u = X, f ^ 



u = X = ±L' ± ^, = 0, aty = ±L', 

satisfies the upper bound 



T < min (^Cia^/V5/3Li/3,C2«V2^^ 



where a = min(ai,a2) and the constants Ci and C2 depend only on K, where K 
is the larger of the supports of $5*^2 divided by a. The constants are independent of 
a, a, L, L' . 
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Proof. We will consider two test solutions, and our upper bound will be the minimum 
of the energies of the two test solutions. 

One test solution is obtained by interpolating between the two boundary conditions 
in the region \y\ < Ka, and connecting this solution to the boundary conditions at 
y = ±L' by a solution with a nearly uniform strain. 

More precisely, we set u = x and 



w~ 



{x,y) = V2, 



aa 



c 



1 X 
2^2L 



+ 



1-C 



1 X 
2^2L 



where ( is as defined in Lemma (5. (i[ Defining 6{x) by 



S{x) = 
we set v{x, y) 



vi{x, y) + V2{x, y) where 



2^ ^1 



\wy{x, ?7)] di] + 2aw~^{x, 0) 



x,y) 



itoo 



[wy{x,r])Ydri + y±-6{x) 



y±ls 



where "i? is a smooth function such that < < 1, {}{x) = 1 for a; < and = 



for X > ^. 



We can estimate the energy of this test solution noting that 6{x) is bounded by 
Ca^a^ by lemma IK3| and using the appropriate rescaling of lemma EIEl to bound the 
energy due to fi. A straightforward calculation gives 



KC 



+ — - — + + 



L3 



+ 



a a L 



This energy corresponds to £i + Em in our heuristic calculation, and this is of course 
reasonable, since we do not have a region corresponding to the self-similar solution 
in this test function. It is also clear that for a <^ L ~ L', the dominant term in 
the energy is cr^a^L/a, and this will be larger than the other terms provided that 
a^a? <^ a'^a^L/a^ i.e in the asymptotic regime a <C a~'^^^a'^^^L^^^. 

Another test solution that we will consider is the following. Let 0^, ip^ be the zero 
asymptotic shift profile with support Kq that we constructed in Lemma [5.41 

We will first assume that ai = 02 = a. Let {u,v^,w^) denote the self-similar 
solution with profile 0^,?/'^, constructed in the same manner as in Theorem. 15. 3[ 
with p(— 1) = p(l) = a. 

We define the test solution by 



u{x,y) 



w 



{x,y) 





— \x\ \ 


(- 


b J 




— \x\\ 


V b / 



u+ 1-C 

W + V2 



L — \x\ 



L 



X 



X 



B 



+ 



1-C 



L + x 
B 
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where ( is as defined in Lemma [5.61 and 6 < L is a length scale we will choose below. 
Defining 6{x) by 



6{x) 



POO 

[wy{x,r])] dr] + [wy{x,r])] dr] + Aw^ {x , 0) 
Jo 



it follows from the construction of w, and of w that 6{x) 
We set f (x, y) = vi{x, y) + V2{x, y) where 

ry 



if L - Ixl > h. 



vf{x,y) 
vtix,y) 



itoo 



[w^{x,r])] d7] + y±-6{x) 



where i? is a smooth function such that < < 1, {}{x) = 1 for a; < |, and = 
for x > ^. Note that v{x,y) = v{x,y) for \y\ < Koa^^^^a^^^L"^^^ and L — \x\ > b. 
We can estimate the energy of this test solution as above and we obtain, 



X < C{K) 



IF 



+ 



^11/3^2^1/3^2/3 



a 
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+ 



a a L 



Remark. The constant C depends on the boundary conditions through the support 
K, but the constant D is independent of K. 

This energy corresponds exactly to the energy 8 in our heuristic calculation. Con- 
sequently, we get 

in the scaling regime a/a ^ 0{1) for the choice h = ot'l^a'^l'^L^I'^ . 



2 2 
a a , 



in the scaling regime aL'^/^a~^/^ ~ 0{1) for the choice b = a e'I^L~'^l^ . D is a 
universal constant, independent of the scaling functions $ and determining the 
boundary conditions. 

□ 



6. Structure of the minimal ridge 

In this section we will derive pointwise and integrated bounds for the ridge width 
p{X) and for the ridge sag W{X,0). Lobkovsky's analysis predicts that these two 
quantities should scale in the same way, and further, the associated length scale is 
not uniform along the ridge j2H] . 

Given a test solution {u,v'^,w'^) that satisfies the boundary conditions, we can 
naturally construct three different x dependent length scales from the solution. They 
are, the inverse curvature 

^0 -no H -1 



ax 



/ _ 2 /"°° 2 

J-oo Jo 
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the ridge sag w'^{x, 0) — w {x, 0), and the support 

k(x) = snp{y > 0\[w^{x, ±y)]'^ + [v^{x, ±y) - yT 5/2]^ + [u{x, ±y)f > 0} 

In our construction for the self-similar test solutions yielding the upper bound, we 
see that a{x),w{x,0) and k{x) all scale in the same way as p{X) — p{x/L). This 
suggests that the structure of the ridge in the ^/-direction is given by a single length 
scale, which depends on x, and further, this length scale is given by our assumed 
scaling for p{X) in the construction for the upper bound. 

Our basic tool will be bounding the stretching energy using only the length scales 
a{x) and w{x, 0) at a given point x. Combining these estimates with our estimates 
for the energy, we obtain pointwise bounds for a{x) and w{x,0). This estimate for 
the stretching energy is obtained in the following lemma. 

Lemma 6.1. Given the profile W{X, .) at a point X e [—1, 1], the stretching energy 
Eg is bounded from below by 



max 



C{l-5) 



VYiw{X,0) + ^ 



y5/2 



6 ^Sp{x)J 



where 



POO rO 

/ W+Y^{X,Y)dY + / WyY^{X,Y)dY 

Jo J-oo 



Proof. The stretching energy is given by 

--nil 

Since the integrand is non-negative, wc have 



W^dX 1 dY 



Eg > max 



y 1 / .1 



WxdX ) dY 



(14) 



W^dXdY 



-Y ^ \J-1 
-Y J -I 

by Jensen's inequality. For each r] e [— F, F], and for all X e (—1, 1), we have the 
elementary inequality 



1 

> max — 
- Y 2Y 



i: 



{W{X,rj)-W{-l,rj)f , {W{X,r])-W{l,rj)) 



l + X 



+ 



> 



(15) 



W^iX,ri) 
2 

W\X,7j) 
1-^2 



1 1 
+ 



l + X 1-X 



1-X 

-2[W'{-l,rj) + W'{l,ri)] 



-2[W'{-l,ri) + W'{l,r})] 
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(±1, 7]) are given by the boundary conditions on the frame. Using the scahng form 
for the boundary conditions, we see that 



(16) 



j W\-l,ri)dri + j W 



All we now need is an estimate for Jf W{^,ri)d^dri. Let p^(X) be given by 



p±(X) = ± 



itoo 



H -1 



W^Y'^{X,Y)dY 



so that [p(X)]-i = [p+(X)]-i + [p-{X)]-\ By Lemma El we obtain 



[ \W^(X,ri)]'^dri > max (1-5) 
Jo z<y,56(o,i)' 



z(w{x,o) + 



P+iX)Z 



+ 



12 



[(3^X)f 



5p+{X) 



Adding the corresponding resuh for W , we get 

f [W{X,r^)fdr^ > max (1-5) zJw{X,0) 

J-Y Z2<Y,5e{0,l) \ 



+ 



P+{X)Z2 



12 



Zo 



+ max (1 — 6) 

Zi<Y,5e{0,l) 



6p+{X) 

( 



+ 



12 



[(3~{X)f- 



Zi 



6p+{X) 



By the matching condition, P^{X) = P~{X) + ^J2. This also imphes that [/5~(X)]^ + 
[/?"'" (X)]^ > 1. Replacing the separate maximization over Zx<Y and Zi < F, by a 
single maximization over Z = Zi = Z2 < F, we obtain 

ljW(X,„)?d„ > ,^m- ,,2(l-.) [z(vr(A-.0) + ^)%|!(2-^)^ 

Combining this with inequalities (fT^ . (fT3j) and (fT^ proves the lemma. □ 
We can now prove pointwise upper bounds for p{X) and W{X, 0). 

Theorem 6.2. {U,V^,W^) is a test solution that satisfies the boundary and the 
matching conditions. Also, I{U,V^,W^) < I. Then, 3Ci, C2, CJ, and C'^ that 
only depend on I such that 

p{X) <Ci(l-X2)2/5 + C[A 
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Also, the ridge sag satisfies the pointwise bound 

W\X, 0) < max (^C[p\X), C',A' (-A_^ , c', (^^^^^) 

Proof. Setting 5 = \,Y = |p(X), Lemma lOl yields 

If p{X) > (2^)1/3 A, it follows that p^{X) - C'A^ > |p^(X). Using this in the above 
inequality, we see that either p{X) < {2C'Y^^A, or 

Combining these two inequalities, we see that 3Ci, C[ only depending on /, such that 

p{X) < Ci{l - X^f/' + C[A. 

□ 

This finishes the first part of the proof. To illustrate the pointwise bounds for 
W'^{X, 0), we begin with a heuristic calculation. If W'^{X, 0) ^ p'^iX), the dominant 
balance in Lemma ffj. II is 

P 

and this gives the characteristic scale Y p. This is the same calculation as above, 
and this gives < < (1 -X^)^/^ as in the previous part. If W'^{X, 0) > p^(X), 
after ignoring the boundary term C'A^F"^/^, the four remaining terms in Lemma lUTTl 
are of orders 

AT- ^ , ^ , y7/2 

iy2(x,0)Vr, l^(X, 0)^3/2, y5/2^ and 

respectively. The dominant balance is between the first and the last terms, and this 
gives the characteristic scale Y ~ [iy^(X, 0)p(X)]^/3. Also, this gives 

^^iy2(x,o)v^~ ^^iy^/=^(x,o)//6(x). 

Rearranging, we get 

We will now make these considerations precise. 
Lemma 6.3. 3^ < oo,C > such that \W{X,Q)\ > Cp{X) implies that 

Vy (w{x, 0) + 1) ' + ^ (i - > c'w'/^x, 0)//3(X), 
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where 



Y 



1/3 



[W\X,Q)p{X)\ 



1/3 



Proof. Let iy(X, 0) = Cip(X), and Y be as defined above. A direct calculation 
shows that 



/ y\2 y5/2 / Y \ 



7/6 



y3 



Cl 



-1/3 



6p(X)J 

1^7/3^1/3 



By taking Ci sufficiently large (> C), we obtain the required inequality. In particular, 
we can take (7 = 8 and C = 1/3. □ 

We can now finish the proof of Theorem 16.21 

Proof. If iy^(X, 0) > C^p(X, 0), combining Lemma f6. II and Lemma [6.31 we obtain 



C 



X2 



If 



it follows that 



[W^\X,^)p{X)) 



I > \ E., > 



1/6 



> 



{W\X,Q)p{X)f\ 



{W\X,Q)p{X)f"' 

C , „s , 1/6 



[W^\X,^)p{X))- 



2(1 -X2) 

Combining all the above considerations, we see that one of the following inequalities 
has to hold 

iy^(x,o) < c[p\x), 

1/4 



'\-x 



2A6\ 1/7 



□ 



This proves the theorem. 

Remark. Our pointwise bounds are not optimal for piX\ In particular, our con- 
struction for the upper bound shows that we can have test solutions with 

p{X) < Cmax {A, e'^/'^il - \X\), (1 - X^^'^) , 

that have a uniformly bounded energy. 

For A <^ (1 — |X|) -C 1, this function is asymptotically (in e) much smaller than 
the upper bound from Theorem 16.21 In the remaining range, i.e, for 1 — |X| < A, of 
for C < 1 — |X|, where C is an 0(1) constant, our pointwise upper bound captures 
the behavior of p{X) in the self-similar construction. 
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Remark. If p{X) did scale like the pointwise upper bound in Theorem 16.21 i.e 
p{X) ~ CiA + C2(l — X^)^/^, then the pointwise bounds for the ridge sag imply that 

\W{X, 0)1 < C[A + C'^il - Xy/'' ~ p{X). 

We cannot obtain lower bounds for p{X) or W{X, 0) purely by energetic arguments 
as the following "pinching" argument shows. For any given point x G {—L,L), we 
can consider the test solution obtained by pinching at x, i.e, we set 

u{x, s) = 0, w{x, s) = V2aa^ ^- j , v{x, s) = a^a"^ ^- j + s, 

where a is a length scale that we are free to choose subject to a > crexp(— e~^/^) 
and ($, \E') is a no-stretch profile with zero asymptotic shift that we constructed in 
lemma 15.41 We use our construction for the upper bound to construct minimal ridge 
solutions in [— L, x] and [x,L]. 

These solutions connect smoothly at x, since our constructions have u{x', s) = 
u{x, s),v{x' , s) = v{x,s) and w{x',s) = w{x,s) for sufficiently small \x' — x\. Also, 
(L + xy + {L — xY < for all p > 0. Combining this with our upper bounds, 

we see that the energy of the "pinched" solution scales in the same manner as the 
upper bound. Since the length scale a can be chosen (essentially) arbitrarily small, 
it follows that the energetics are not enough to give us a pointwise lower bound on 
the ridge-sag W{X,0) or the the ridge-width p{X). 

We need these pointwise lower bounds to prove rigorous scaling results for p{X) 
and W{X,Q).To obtain such results, we need to use the fact that the solution of 
interest is a minimizer, i. e, the first variations vanish. This type of analysis is carried 
out for the structure of an Austenite-Martensite boundary in ITU] . A similar 
analysis is possible for the minimal ridge, and we will present the details elsewhere. 

7. Discussion 

We conclude our discussion by indicating some of the issues/open problems relating 
to thin elastic sheets in general and to minimal ridges in particular, and in the process 
we point out the relevance of our results to some of these questions. 

We have proved rigorous scaling laws for the energy of a single minimal ridge with 
a geometrically linear FvK ansatz. A natural question is the extension of these results 
to the Nonlinear FvK energy, and also to full three dimensional elasticity, as in 
It is easy to extend this to a mixed energy functional where the bending energy is 
treated in a geometrically linear fashion, but we use the "full" energy W2D for the in- 
plane stretching. Extending this analysis to the Geometrically nonlinear functional, 
or to full three dimensional elasticity will require some new techniques |4j. 

Another problem is to show that the scaling laws also hold for a real crumpled 
sheet, where the forcing is not through clamping the boundaries to a frame, but rather 
through the confinement in a small volume. In this case, there are interesting global 
geometric and topological considerations, some of which are explored in Refs. 
and |16j. As Lobkovsky and Witten [29j argue, the boundary condition that the 
deformation goes to zero far away from the ridge implies that the ridges do not 
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interact with each other significantly. The ridges can be considered the elementary 
excitations of a crumpled sheet. 

More precisely, we have constructed ridge solutions with zero asymptotic shift, 
that are exactly strain free on the boundaries. These solutions give "non-interacting" 
ridges and patching these solutions together, it is possible to construct test solutions 
for a sheet confined inside a sphere. This gives us upper bounds which scale in the 
same way as the energy of a single ridge. 

To show the corresponding lower bound, we have to show that confinement actually 
leads to the formation of ridges, and that the competition between the bending and 
the stretching energy for this situation has the same form as in lemma In this 
context, we expect that global topological considerations, as well as the non self- 
intersection of the sheet will play a key role in the analysis, as they do in the analysis 
of elastic rods (one dimensional objects) pUj. 

As we note above, the blistering problem is described by the same elastic energy 
(Eq. Q), but with different boundary conditions. Our results show an interesting 
contrast with results for the blistering problem. Ben Belgacem et al. have shown that 
^ , for an isotropically compressed thin film, the energy of the minimizer satisfies 

cA^/VL < J < ca^/Vl, 

where L is a typical length scale of the domain, and A is the compression factor. A 
construction for the upper bound strongly suggests that the minimizers develop an 
infinitely branched network with oscillations on increasingly finer scales as a — 0. In 
contrast, our results indicate that the energy of a minimal ridge satisfies 

^^5/3^1/3 < J < Ca^/s^Vs^ 

and the energy concentrates in a region of width a^^^L'^^^. This shows that the nature 
of the solution of the variational problem for the elastic energy in Q depends very 
strongly on the boundary conditions. In particular the very nature of the energy 
minimizers is different for the two problems - For the blistering problems, as a — 
the minimizers develop a branched network of folds refining towards the boundary. 

Finally,one would like to prove F-convergence and find the F-limit for the 

elastic energy as cr 0. The difference in the scaling of the energy minimum for 
the minimal ridge, and the blistering problem shows that the F-limit of the elastic 
energy depends crucially on the imposed boundary conditions. 

The analysis in this paper only pertains to situations where the configuration of 
the sheet is either smooth, or consists of a finite number of minimal ridges. More 
precisely, the sheet configurations : 5 — > is piecewise smooth, strain- free a.e., 
has gradient D(j) in BV, and the singular support of D'^(j) lives on a finite union of 
straight line segments. For the boundary conditions that admit such configurations, 
our analysis suggests that the asymptotic energy is on the scale a^^^ as a 0. 
Further, if the following limit exists, the is necessarily follows that 

T[0] = F-lhna-^/3j = C$:af/f, 
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where Ij is the length of the jth segment in the singular support of -D^0, and aj is 

1/3 

the jump in D(j) across the segment aj. Note that, because of the Ij dependence, 
the F-limit cannot be written as the integral with respect to the Ti^ measure on the 
singular support of D'^cj), of a local energy density, which only depends on aj. 

The reason we get a dependence instead of a linear dependence in Ij is that the 
ridges have a nonuniform structure along their length for any a > 0. Any formulation 
of the F-limit should therefore incorporate a "hidden" variable which reflects the 
nonuniform structure, although, this non- uniformity is no longer detectable in the 
limiting a —>■ configuration. A natural candidate for this variable is the scaled 
ridge width p{x) (= the inverse curvature a{x) defined as in Sec. where x is a 
coordinate along a ridge. We expect that the F-limit can be written as an integral 
of an energy density with respect to measure on the defect set, where the energy 
density depends on a and also on p(x) and derivatives p'{x) and p"{x). In fact, 
Eq. fjl2|) strongly suggest that the F-limit for the energy of a single ridge can be 
written as 



= cy «;/^nf 

^ ' Pi-) Jo 



{[p"ix)pix)r + [p'ix)]') p{x) + ^ 



dx 



p{x)_ 

where x is a coordinate along the ridge, and the infimum is over all smooth functions 
p that vanish at both the endpoints of the ridge. 

We hope this paper spurs further investigation of the question of the F-limit of 
the elastic energy. This question is very much open, there are no proofs for either 
F-convergence, or of our conjectured structure for the F-limit. 
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